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ABSTRACT 

A charged particle passing through or by an optical fiber induces emission of 
light guided by the fiber. The formula giving the spontaneous emission ampli- 
tude are given in the general case when the particle trajectory is not parallel to 
the fiber axis. At small angle, the photon yield grows like the inverse power of 
the angle and in the parallel limiting case the fiber Cherenkov effect studied by 
Bogdankevich and Bolotovskii is recovered. Possible application to beam diag- 
nostics are discussed, as well as resonance effects when the particle trajectory 
or the fiber is bent periodically. 

keywords optical fibers, monomode fibers, cherenkov radiation, transition ra- 
diation, beam diagnostics 

1 Introduction 

A charged particle passing through an optical fiber of large enough (a ^S> A) ra- 
dius produces Cherenkov radiation (CR) inside the fiber if its velocity is above 
the Cherenkov threshold f / %fe. Photons which are emitted at sufficiently small 
angle (cos 9 > 1/ ^/e) with the fiber axis undergo total reflection and can be 
guided by the fiber toward a photon detector, making a position sensitive par- 
ticle detector PQ. 

Transition Radiation (TR), emitted when the particle crosses the fiber sur- 
face, is another mechanism producing guided light in the fiber. Such a mecha- 
nism has been suggested to produce X-rays in capillaries ( "channeled transition 
radiation" j2j). It has no threshold in velocity. Above the Cherenkov threshold, 
adding naively Cherenkov and Transition Radiations makes double counting. 
Indeed, Transition Radiation from the entrance and exit faces of a dielectric 
slab, added coherently, contains Cherenkov radiation as a singular part. Only 
for broad enough fibers can one apply the Cherenkov formula, neglecting the 
edge effects. 

The Transition Radiation formula, derived for the ideal case of an infinite 
flat boundary, can be applied to the case of a lamellar optical guide. However it 
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is inadequate for narrow (a ~ A) cylindrical fibers which have only one or a few 
number of modes, because of the strong curvature of the surface. It is precisely 
this case (a ~ A) that we will study here. 






Figure 1: Standard mechanisms of light emission. Top: Cherenkov radiation. 
Bottom: Transition radiation 

Photon emission inside the fiber is also like the cavity radiation considered 
by Nitta and Miyazaki jjJJ. They consider rectangular cavities of finite size in 
all dimensions, where the frequency spectrum is discrete. In our case we have 
one quasi-infinite dimension and continuous spectrum. 

If the particle passes close to the fiber without touching it, its Coulomb field 
can be diffracted, giving rise to Diffraction Radiation. This radiation is mainly 
outside, i.e. in the "free" modes of the fiber, but nothing forbids some part to 
be kept "inside" . This phenomenon will also be treated on the same footing as 
the case where the particle crosses the fiber. 

2 Basic radiation formula 

We take the fiber axis as z axis. The quantized electromagnetic field A op in 
absence of the particle but in presence of the fiber can be expanded in proper 
modes of the longitudinal momentum k z (continuous), total (orbital + spin) 
angular momentum M = L z + S z (discrete) about the z-axis and radial quantum 
number v (continuous or discrete): 
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The radiation gauge ^4° = has been chosen. We have quantized the field in 
an infinitely long cylindrical box, so that k z is continuous and v is discrete not 
only for guided modes, but also for free ones. The photon energy to depends on 
k z , M and v through a dispersion formula 



uj = uj(M, k g ,v). 



(1) 



Fig. 2 shows the phase velocity — o-i/k z of the lowest mode (M — ±1, v = 1) 
called if -En and the external fraction of the mode power, as a function of u. 




at . a 



Figure 2: phase velocity — uj/k z and external fraction of the power for the 
HE\\ mode. 

The time-independent annihilation and creation operators obey the commu- 
tation rules 

= 2ir5(k z - k' z ) 5mm- <W- (2) 

A particle of charge e is crossing or passing by the fiber along the (straight 
or curved) trajectoty X(t). The amplitude R of spontaneous emission in the 
mode k z , M, v is given by 

pt — OO 

R*m^u = e dX(t). [A M ,k,Az,v) e izK -^\ (3) 

J t — — OQ 

= — [ dX ■ E M ,k,Ax, V) e izk ^ (4) 
iuj J 



0,M,k z ,v, a M',k' z y 
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where E = iuiA is the associated electric field, according to E = — dt A — VAo 
The resulting photon distribution is 



dN M ,k z ,u = -77— \B,M.kz,v\ 2 ( 5 ) 

3 The field modes 

The fiber we consider consists in a dielectric rod of homogeous index surrounded 
by vacuum (no clading). A guided mode is such that 

to < k z < <jj\fe = u>n. 

Outside the fiber, it has an evanescent wave ~ e~ Kr corresponding to an imagi- 
nary transverse momentum in. Inside the fiber, the photon has real transverse 
momentum q. q and k are given by 

q(M,k z ,v) = {ecu 2 - kl} 1 ' 2 K(M,k z ,v) = {k 2 z -cu 2 } 1/2 . 

The longitudinal components E z and B z of the electric and magnetic fields 
have longitudinal spin S z — so that their orbital angular momentum L z = M. 
Using cylindrical coordinates (z, r, (f>) we can therefore write 

Ei M ' k ""\x,y) = i e iM * Mr), B^ v \x,y) = e iM * h,(r), 

Inside each homegeneous region (i.e., core and vacuum) f z and h z obey the 
same differential equation 

r 1 M 2 

d r+- d r-^T+ k T(r) 

where 

k T (r) = w 2 e(r) — k 2 z — q 2 inside, — k 2 outside, 

e(r) = e inside, 1 outside. 

The solutions of iJBJ are the Bessel functions Jm and Km- We have the limiting 
conditions: 

• f z and h z are regular at r = 0. Therefore f z (xh z cx Jh[{qr) for r < a, 

• f z and h z are continuous at r = a, 

• f z and h z vanish at r — > 00. Therefore f z och z oc Km{kt) for r < fl- 
it follows that f z and h z are proportional to each other (it would not be the 
case if there were more than one boundary) and we can write: 

fz(r) = c E i/>(r), h z (r) = c B i>{r), (7) 



(6) 
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ip(r) = Jm(qt) inside, = ck Km(kt) outside, ck — 



•hiiqa) 



K M (na) 

ip is continuous at r — a, but not -0'. So 10 cannot be treated as an ordinary 
Schr0dinger equation. 

The transverse components and Bx can be expressed in terms of the 
radial and azimutal vectors, e r = cos <fi x + sin <fi y and = cos 4> y — sin cj> x, 
times L z = M radial functions: 

Er = e m + (f r (r) e r + U(r) e*) 

B T = e lMfl > (h r (r) e r + h^r) e*) 

They can also be expressed in terms of the S z = ±1 eigenvectors — (x±iy)/2, 
times L z — M q= 1 radial functions: 



Et = e i[M-l)4 f {r) e+ + e i{M+l)t /+(r) e - 

^B T = e < M ~ x ^ h^{r) e+ + e l < M+1 ^ h+(r) e 



(8) 
(9) 



with 

f± = f r ± if<f> , -ih± = h r ± ih 4 
The Maxwell equations lead to (see Appendix): 



f± = 



-k z c E ±ui c b 



-fc z c_b ± ^e(r) c E 
4(r) 

Owing to a property of the Bessel functions, 



r 



r 



Jw(a;) ± —J M (x) = ±J Mt i( x ) 
x 

one obtains the final expressions 

f±(r) = cfj J M±1 (qr) (r < a), 
h±{r) = c±j J M ±i{qr) (r < a), 



M 

K' M (x) ± — K M {x) = -K Mt i(x), 
x 



cf K K M ±i(nr) (r > a), 
c hK k m±i(kt) (r > a), 



(10) 
(11) 



with 



-hj 



-k z c e ±uj c B 

q 

-k z c B ± toe c E 

q 



-k z c E ±u! c E 



-fK 



J hK 



-k z c b ±oj c E 



c K - 



The continuity of h z , h,p, f z , and e(r)/ r at the fiber surface leads to 



CE 



-MQ 



J' M {u) , K' M (w) 



uJm(u) wK m {w) 



MQ 



eJ' M {u) , K' M (w) 



uJm{u) wK m {w) 
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where u = qa, w = na and 

Q 

from the two expressions of cb/ce we obtain 



1 


1 " 




e 


1 " 




u 2 


v — 

ur 


UJ 









J' M (u) , K' M (w) 



lJ m {u) wK m {w) 



sJ' M ( u ) , k 'm( w ) 



uJ m {u) wK m (w) 



M 



This equation, together with 

u 2 + w 2 — (e - 1) uj 2 a 2 and m 2 + £w 2 



(e-1) k 2 z a 2 



1 

(12) 



(13) 



determines the dispersion relation Q between w and 

For optical hbers used in telecommunications, e = e\jei is the relative per- 
mitivity between the core and the clad, and is close to unity. Then E x and E y 
can be treated like independent scalar fields and e is replaced by 1 in l|12|). In 
our case, we have no clad and e is large, for instance s = (1.41) 2 ~ 2 for fused 
silica. So we cannot make such an approximation. 

4 Normalization of the mode 

The commutation relations © imply that the mode A is normalized to unit 
linear photon density, i.e., a linear energy density 

dW SM 

— {A} =w. 

dz 

This fixes the absolute value of the coefficients Ce and cb- However the elec- 
tromagnetic energy density has not a simple expression if there the medium is 
dispersive. On the other hand the power is given by the simple expression 



P = ^{A} = J dxdy 5R{E* x B} 2 



where 5i {E* x B} is the Poynting vector. In terms of the functions introduced 



P=l j2nrdr^{r_(r) h_(r) - f* + (r)h+(r)} 
= 7T J r dr {cjj clj J 2 m _ x (qr) - cjj c+j (gr)} 
-7T J r dr \c~ fK cy iK Kl^^nr) - cj K c\ K K 2 M+l {nr)} 



(14) 



(15) 



Using the two expressions of the group velocity 

P dijj 



dW/dz dk z 
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we can transform (J5J as follows 

du 



dN 



M,k z ,v 



2ttv 



. .2 U) dLO 2 

\RM,k.,v\ = 7; — 77 -Rm.L.i/ • 
Z7T_T 



(16) 



This formula is invariant under an overall rescaling of the fields ce and cb since 
both P and |i?| 2 are quadratic in the fields. 



5 The mode emission amplitude 

We assume that the particle follows the straight trajectory 

X = b + tV, b= (6,0,0), V=(0,V T ,V z ) 
The spontaneous emission {amplitude}* becomes 



rL M,k z ,v 



e 

2uj 



e 



dy 



dy 



f-(r)-e** / + (r)+2-^(r) 



V, 



v T - 



One arrives at the real expression 



RM.k z ,u = - dy{ cos[sy + (M - 1)0] /_(r) - cos[sy + (M + 1)0] /+(r) 



+2(I4/Vt) cos( S2 / + M0) / z (r) } (17) 



where r = Vp + j/ 2 , = tan ^y/fe), s = (V z k z - u)/Vr', /+, /- and f z are 
given in (|7ITrj|) . 

For small transverse velocity, the third term of the integrand becomes large, 
although the longitudinal component f z is generally smaller than the transverse 
ones, and the photon spectrum can reach high values. On the other hand the 
amplitude is expected to decrease very fast when the dimensionless parameter 
sa gets much larger than unity, because the integrand oscillates many times 
in the important field region \y\ < a. From these considerations, at small Vr 
one expects an almost monochromatic peak of height cx (1/Vr) 2 centered at 
frequencies such that s = 0, i.e., the electron and the wave travel at the same 
longitudinal velocity: 

tu/k z = v^ = V z (18) 

("fiber Cherenkov condition"). The width Aw of this peak is given by the 
condition a As ~ 1, where 



a As = \Auj - KAL 



a 



|1 - Vi/v g \ 



a Alo dv^/vcf) a Aw 



V T 



dw/uj Vt 



1 (19) 



Given v z , the Cherenkov frequency lock can be obtained from the upper curve 
of Fig. 2. Also the inverse of dv$j{v§ ckv) is the interval between u>ch and the 
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intersection of the tangent to this curve and the horizontal axis. Aw being 
proportional to Vt, after integration over to one obtains a net photon number 
oc 1/Vt- At small Vr, one approaches the case when the electron runs parallel 
to the fiber over a long distance L 0], with L ~ a/Vr- The spectrum becomes 
discrete, given by the "fiber Cherenkov condition" (JTSJ. 



6 Results and discussion 




Figure 3: Photon yields in the HEu mode for four types of the particle trajec- 
tory. 

We have considered fibers of fused silica with n=1.41. The dimensionless 
photon yields LodN/duj is plotted in Fig. 3 for two impact parameters, b — 
0.2 a (penetrating trajectory) and b = a (tangent trajectories), and two particle 
velocity vectors, (v z ,vt) — (0.88,0.1) (large angle) and (v z ,vt) = (0.85,0.5) 
(small angle) . The angular momentum M of the mode has been chosen to have 
the same sign as the particle's one. This is expected to be more favorable for 
radiation than the opposite sign case. 

The spectrum from the penetrating trajectories are harder than from the 
tangent (or fully external) trajectories. This can be due to (i) the discontinuity 
of the fields at the fiber surface, (ii) the fact that the evanescent field is less 
important at high frequency. 

In the large angle - penetrating case, the yield is of the order of a = 
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e 2 /(47r) = 1/137, as more or less expected. For the large angle - tangent case 
it is two orders of magnitude smaller. Note the peak at a relatively small fre- 
quency, where the wave travels mainly outside the fiber (see Fig. 2). At still 
smaller frequency, the wave becomes too much diluted outside the fiber, which 
explains the vanishing yields in the four cases. 

Much bigger yields are obtained in the small angle cases, due to the increas- 
ing contribution of the longitudinal field discussed in the last section. Ea. (|18fl 
should predict a peak at ua ~ 2. In the tangent case, this peak is shifted to 
the left, probably for the reasons (i) and (ii) said before. In the penetrating 
case, instead of a peak, we have a dip at this very place. This is a peculiarity 
of the odd M modes when b is small. If b = 0, then </> in (|17fl is either —tt/2 or 
+7r/2 and, at the Cherenkov point (s = 0), cos(sy + Mcf>) is zero in the whole 
integration range. 

7 Outlooks 

The production of coherent light by charged particle in a fiber may have appli- 
cations in beam size measurement. Scintillating fibers yield more photons, but 
may be less hard to radiation than, e.g. fused silica. One can also compare with 
optical transition radiation and diffraction radiation devices. In our case, the 
synchrotron light from an upstream magnet or the diffraction radiation from a 
collimator does not make any background. It is just scattered, but not captured 
by the fiber, due to the orthogonality between the free and guided modes. 

Unlike for optical transition radiation and diffraction radiation we do not 
expect a logarithmic increase at large Lorentz factor 7. Such an increase is usu- 
ally due to the expansion of the photon cloud accompagnying the particle, but 
such photons are quasi-real, therefore almost orthogonal to the guided modes. 

Eq.QJ can be applied to cases of bent particle trajectory or, using the adi- 
abatic approximation for the modes, of bent fiber. Successive encounters will 
give rise to interference phenomena. A straight fiber inside an electron undu- 
lator will get filled by monochromatic radiation, to and k z being related by the 
resonance condition 

(k z -uj/(V z )) A = k z A-uoL/V = 2mn (20) 

where A is the undulator period and L = AV/ (V z ) the lenght of one period of 
trajectory. Similarly, a particle in straight motion close to a periodically bent 
fiber will induce monochromatic light in the fiber. This is a kind of internal 
Smith-Purcell radiation. If the fiber draws a planar periodic curve, the resonance 
condition is 

k L L - ljA/V = 2rmr (21) 

where now A is the spatial period of the fiber, L its corresponding arc length 
and the wave number along the fiber. If the fiber draws an helical curve, the 
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phase velocities of left- and right-handed polarisations in the fiber are split and 
the condition becomes 

k L L - ujA/V = 2mir ± (j) B (22) 

where <pg, called Berry phase, is equal to the solid angle of the cone drawn by 
the tangent to the fiber [Hj. 
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